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Entanglement between large numbers of quantum modes is the quintessential resource for 
future technologies such as the quantum internet. Conventionally, the generation of multimode 
entanglement in optics requires complex layouts of beamsplitters and phase shifters in order 
to transform the input modes into entangled modes. Here we report the highly versatile and 
efficient generation of various multimode entangled states with the ability to switch between 
different linear optics networks in real time. By defining our modes to be combinations of 
different spatial regions of one beam, we may use just one pair of multi-pixel detectors in order 
to measure multiple entangled modes. We programme virtual networks that are fully equivalent 
to the physical linear optics networks they are emulating. We present results for N = 2 up to 
N = 8 entangled modes here, including N = 2, 3, 4 cluster states. Our approach introduces the 
highly sought after attributes of flexibility and scalability to multimode entanglement. 
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Multi-partite entanglement is not only of fundamental 
scientific interest, it is also the key ingredient for quantum 
information technologies^"^. In optics, several impressive 
demonstrations of multi-partite entanglement have been shown 
recently including an 8-photon cluster state^ and a 9-mode state 
used for error correction^. However, these schemes tend to employ 
one detection system per entangled mode/qubit, which intro- 
duces a lack of flexibility and is detrimental to its scalability. These 
optical setups are built to produce one set of outputs or to per- 
form one given protocol; in order to change the output the optical 
hardware itself must be modified. 

Currently, the well-established recipe for generating entangle- 
ment using continuous wave laser beams is to mix squeezed modes 
of light together at beamsplitters. It is possible to create iV-mode 
entanglement given a network of iV- 1 beamsplitters with N input 
modes, even with less than N squeezed modes^. 

In our scheme, we copropagate all possible spatial modes of light 
within one beam. Entanglement between copropagating modes in 
one beam has been previously demonstrated with spatial modes^'^, 
and also in the frequency domain^^. In the current work, we 
radically expand the idea of one-beam entanglement by introduc- 
ing the notion of emulating linear optics networks, by programming 
virtual networks that mix together different spatial regions of the 
light beam. These software-based networks calculate the precise 
weighted combinations of the spatial regions required to emulate 
the physical networks. This is possible because the linear optical 
components in a typical network simply perform reversible opera- 
tions, and can be represented by unitary matrices. It is worth stat- 
ing explicitly that the entangled spatial modes that we produce are 
event-ready and unconditional before the detection process. The 
real-time virtual networks allow us to match the detection basis 
to the desired spatial mode basis contained within the beam, analo- 
gous to shaping a reference local oscillator beam. We report here 
the ability to switch in real time between desired output states in 
a one-beam entanglement system using just one detection scheme. 
We programme virtual networks for 2, 3, 4, 5, 6, 7 and 8-mode 
entangled states. As a further demonstration of the versatility of our 
setup, we produce linear 2, 3 and 4-mode cluster states, which are 
highly entangled graph states garnering attention for their potential 
in quantum computing^. The cluster state measurements demand 
more stringent squeezing requirements than non- cluster insepara- 
ble states. 

Results 

Measuring spatial modes. By employing custom-made multi- 
photodiode-homodyne-detectors (MPHD) that each contain an 
array of eight photodiodes (see Fig. 1), we detect the light in eight 
spatial regions and assign individual electronic gains to each one. 
The linear combination of the eight gain-adjusted photocurrents 
constitutes the measurement of one mode. 

More generally, we can express the measurement process of a 
complete set of spatial modes in one beam by the following: 

a = Ui (1) 
= ^net^ini (2) 

where si = (ai,-'-,a^) is the set of N measured modes projected 
by the Nx8 unitary matrix U acting on the eight homodyne- 
subtracted photocurrent operators i = {iiy---yi^)^ . uf^ is an JVxS 
matrix made up of the top N rows of U^^, the orthogonal 8x8 uni- 
tary matrix that recovers the important set of eight unmixed spatial 
modes that span the input basis (see Methods). Input modes are 
then mixed via (7net» which emulate linear optics networks, given by 
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Figure 1 1 Experimental setup. Squeezed light is prepared and combined 
in 'Squeezers' with a piezoelectric transducer (PZT) controlling the phase 
between the two squeezed modes, locked in quadrature. Vacuum modes 
(vac) copropagate so that the beam exiting 'Squeezers' and entering 
'Detection' contains eight measurable spatial modes. MPHD is used to 
measure the quadrature amplitudes of the beam in eight different regions, 
in 'Detection'. Local oscillator (LO) gives a reference to phase quadratures. 
A PC is used to calculate electronic gain functions via the notion of 
'Virtual networks'. The detected beam is then projected onto a basis of 
measured modes (see equation (1)). Inset: FM generation; half of the wave 
is phase retarded by half a wavelength, flipping the electric field amplitude. 



the NxN matrix: 
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where g R. 

This allows us to uniquely define a mode by the eight real 
numbers in the nth row of 11, which we will label as the mode's gain 
vector G„, such that = G„i. Therefore each spatial mode we meas- 
ure, whether belonging to the input basis or an entangled mode 
basis, is defined by a unique pattern within the light beam. These 
spatial mode patterns, represented by Gaussian profiles modulated 
by respective electronic gains G„, are shown visually in Fig. 2, while 
the detection stage of Fig. 1 shows how we implement this experi- 
mentally. The spatial modes are orthogonal to each other, span- 
ning a basis so that the independent measurement of each mode is 
possible^ 

Input basis. We create two amplitude squeezed modes using optical 
parametric amplifiers (OPAs). The first mode is converted to a flip 
mode (FM) by phase delaying half its beam by half a wavelength, k 
(see inset of Fig. 1). The FM is overlapped in quadrature with the 
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Figure 2 | Spatial mode patterns. Measured modes are defined by spatial patterns of electric field amplitudes. Shown in the top right box is an example of 
how the spatial mode pattern for 05 is matched by applying eight electronic gain values (G5) to the detected Gaussian profile (i). The basis of input modes 
ay..a^ is shown in the middle row (see Methods). The arrows represent a mapping via the virtual networks L/^gt (blue) and L/^gt (violet) onto the respective 
bases of entangled modes; the top row shows the symmetric EPR or 2-mode basis, while the bottom row shows the 8-mode basis. Superscripts denote 
mode basis and subscripts denote mode number. The spatial mode bases for N = 3 to N = 7 are not shown for brevity. 




Figure 3 | Multimode entanglement via emulated linear optics networks. Squeezed light and vacua are mixed together using unitary operations in 
order to produce entangled mode states. Unless otherwise stated beamsplitters are 50% reflective. Other reflectivities are denoted by splitting ratios; 2:1 
denotes a 33% reflectivity, and 3:1 denotes a 25% reflectivity. Superscripts denote mode basis and subscripts denote mode number; there is a one-to- 
one correspondence between spatial mode bases shown here and those shown in Fig. 2 ({G/},{a,^}, {af }). (a) The emulated linear optics network used to 
measure 2-mode EPR entanglement (^net)- 8-mode entanglement via a calculated concatenation of beamsplitter and n phase shift operations (^net)- 
The dots between a and b imply virtual networks for N = 3 ... 7, again not shown for brevity. 



Gaussian mode (GM) output of the second OPA upon reflection 
of its output coupler^^. These two squeezed modes are the first two 
modes of what we refer to as the input basis; and Six copropa- 
gating vacua modes are measured by calculating G„ vectors that are 
orthogonal to both and These vacua modes (labelled ao^...a<^ 
complete the input mode basis (see middle row of Fig. 2). Meas- 
uring these modes amounts to matching the detection basis by 
following equation (1) and setting (7net-^ (see Methods). 

Entangled mode bases. Programming a virtual network amounts 
to calculating the precise expression for (7net- The unitaries we have 
access to in programming the virtual networks are beamsplitters 
and TT phase shifts. (7net is the concatenation of all of these unitaries 
that make up a linear optics network. The tt phase shift is equivalent 
to multiplying a by - 1. Note that arbitrary phase shifts are forbid- 
den as each measurement naturally corresponds to detection at a 
fixed phase defined by a shared reference beam, the local oscillator. 
Importantly, optimal virtual networks are calculated allowing for 



optimization of beamsplitters owing to asymmetries in the squeezing 
levels of input modes. 

The most intuitive virtual network we create is the 2-mode Einstein- 
Podolsky-Rosen (EPR) state^^ shown in Fig. 3a with two squeezed 
inputs. Here, we engineer spatial mode patterns that have no spatial 
overlap; the left half of the beam is entangled with the right half (see 
the iV=2 panel of Fig. 2). Entangled modes belonging to other bases 
share spatial overlap but are nevertheless spatially orthogonal. 

Spatial modes measured in an entangled mode basis are given a 
superscript iV to distinguish them from modes in the input basis; 

and ^2 represent the two modes spanning the iV= 2-mode EPR 
basis. See the Methods section for details on how we create virtual 
networks for each of the iV=2, 3, 4, 5, 6, 7 and 8-mode bases. In 
general, we construct networks pertaining to iV modes by concate- 
nating iV- 1 virtual beamsplitters with vacua on unused input ports 
(see Fig. 3b). This is a highly efficient approach to create multimode 
entanglement as the arduous tasks of mode matching and alignment 
are replaced with the ease of programming. 
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Figure 4 | Noise variance measurements of the spatial modes, (a) x quadrature measurements of the input mode basis. The squeezed ([Ax]]^) is shown 
in red and antisqueezed ([Ax2]^> is shown in blue. The x quadrature variances of the six vacua modes are measured to equal quantum noise (shown by the 
six coloured traces at OdB). (b) p quadrature measurements. The antisqueezed (lApi']^) is shown in red and squeezed ([Ap2]^> is shown in blue. The p 
quadrature variances of the six modes are again measured to equal quantum noise, confirming they are vacua (shown by the six coloured traces at OdB). 
(c) These variances show the x quadrature correlations between modes as in the first half of the left hand side of equation (4). Every column shows N-^ traces of 
X quadrature correlations below shot noise, as well as the blue shot noise trace (OdB) normalized to two units of vacua. Each green trace shows ([A(xf -X2 )^ > 
for each N-mode basis. Each new colour represents the other N-1 variance correlation traces of equation (4). (d) Correlations between measured modes 
in p quadrature, second half of the left hand side of equation (4). Each green trace now shows ([A(pf +P2 +g3P3 +... + g/vP/!/)]^. The traces overlapping 
show that each pair of modes is entangled with the same strength as any other pair of modes, a result of optimizing for symmetry in the virtual networks. 



Cluster states. Attracting attention for their potential in one-way 
quantum computing schemes, cluster states are a type of highly 
entangled Gaussian graph state They satisfy the quadrature 
relation (pa-^beN ^^)^0- As infinite squeezing would require 
infinite energy and are thus unrealizable, one is limited to the pro- 
duction of approximate cluster states in the laboratory, and there 
have been demonstrations of up to 4-mode continuous -variable 
cluster states thus far^^'^^. 

In order to measure cluster states in one beam, we must be able 
to access the correct quadratures of each entangled mode. Here we 
measure 2, 3 and 4-mode linear cluster states; however, measuring 
arbitrary cluster shapes would require modifying the optical setup 
(see Methods). 

Discussion 

Each spatial mode is characterized by the continuous -variable quad- 
rature operators x and p of the electric field operator. The x and p 
variance measurements of the eight modes in the input basis are 
shown in Fig. 4a,b. Here, {[Axq^]^) = {[/^Xi]^) = - 4.3 ± 0.05 dB and 
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([ApFM]^) = ([Ap2]^)=-3.7±0.05dB below the quantum noise limit, 
and the variances of the vacua are verified to equal quantum noise. 

In order to verify entanglement between measured modes, 
we use the well-established van Loock-Furusawa inseparability 
criteria^^. For an iV-mode entangled state, it is sufl&cient to satisfy 
N-1 inseparability inequalities: 

I: 

([A(Xi-X2)f) 

+([ A(Pi + P2 + g3p3 + - + gNpN )f)<h 
N-l: 

{[A{x^_i-x^)f) 

HlMgiPi + ... + gN-2pN-2 + Pn-1 + Pat ) ]^ ) < 1 . 
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Table 1 1 inseparability of entangled modes. 
N 1 II III 


IV 
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VI 


VII 


Average 


2 0.39 














0.39 


3 0.56 


0.56 












0.56 


4 0.64 


0.63 


0.64 










0.64 


5 0.69 


0.69 


0.70 


0.70 








0.69 


6 0.73 


0.73 


0.75 


0.74 


0.74 






0.74 


7 0.77 


0.78 


0.77 


0.76 


0.77 


0.77 




0.77 


8 0.79 


0.79 


0.78 


0.81 


0.79 


0.80 


0.79 


0.79 


Data are based on the van Loock-Furusawa criteria, 
This verifies entanglement of the N modes. 
Uncertainty is +0.01 in all cases. 


given in equation (4). Each row shows that for a basis of N quantum modes, the N- 


-1 values obtained from 


quadrature variances are well below 1. 



Table 2 | inseparability of cluster states. 
N 1 


II 


III 


IV 


Average 


2 0.39 








0.39 


3 0.49 


0.70 






0.59 


4 0.79 


0.67 


0.84 




0.76 


5 0.79 


0.67 


1.10 


1.18 


0.93 


For each /V-mode cluster state, the N-1 inseparability conditions given in equations (8)-(ll) must be satisfied. The partial violation for N=5 means that the state is not fully inseparable, and hence not 
a cluster state. 

Uncertainty is ±0.01 in all cases. 



with free parameters gf, to be optimized for maximum inseparability. 
We have omitted the superscript JV here for clarity, as the above 
holds for any mode basis. The subscripts n of x and p here indicate 
the nth mode in the iV-mode basis. Table 1 summarizes the meas- 
ured degrees of inseparability for all N-1 inequalities in each 
iV-mode basis, given in roman numerals. Table 2 summarizes the 
more stringent inseparability required for unweighted cluster states 
(all homodyne gains {gj} are set to 1). The relevant inseparability 
inequalities satisfy the cluster state quadrature relationship writ- 
ten in the form of equation (4), and are written out explicitly in the 
Methods for reference. 

The terms in equation (4) measure the degree of correlations 
between any two modes in a given basis. For the modes to be 
inseparable, each of these correlation variances (correlations in the 
X quadrature and anticorrelations in the p quadrature) must be in 
the quantum regime, that is below the normalized quantum noise 
of two units of vacua. Figure 4c,d shows this to be the case in our 
experimental measurements. Although we are limited here to eight 
modes owing to our detector, this scheme is scalable to higher num- 
bers of mode entanglement even without increasing the number 
of squeezing resources, as shown in the simulation traces of Fig. 5. 
As we increase the simulated number of modes in the basis up to 
30, the degree of inseparability approaches the classical bound of 1 
owing to the vacuum noise penalty for each additional unsqueezed 
mode input. Entanglement is shown to hold here, however, even 
with current squeezing levels. Importantly, there is no loss incurred 
during the transformation of the squeezed input modes into a set 
of entangled modes, as can be seen by the agreement of the theo- 
retical predictions and the experimental values of Fig. 5a. This 
equates to perfect mode matching at every virtual beamsplitter. 
Figure 5b explores how inseparability scales with different squeezing 
levels. Measuring a larger number of inseparable modes experimen- 
tally requires only an increase in the number of photodiodes in the 
MPHD, and importantly no modification of the optical setup. Note 
that this is not true for cluster states, and the number of squeezed 
inputs must be increased accordingly. 

For the special case of N=2, optimal EPR entanglement^^ 
is measured to be 0.58 ±0.01. Optimizing for the beam-splitter 



reflectivity^ we find that owing to the slight asymmetry between 
input squeezing levels, the optimal beam-splitter ratio here is not 
50%, but rather 48.8%, leading to a very slight improvement over the 
symmetric network. Each unique beam-splitter reflectivity changes 
the mapping of (7net such that formally the beam of light contains an 
infinite number of mode bases. The versatility of our scheme comes 
from being able to match the detection basis to a network that has 
been optimized for an arbitrary set of inputs. 

The entanglement demonstrated in the current work allows for 
such protocols as quantum teleportation^'^^'^^. To perform complex 
protocols such as one-way measurement-based quantum compu- 
tations^'^^'^^, we need to increase the degrees of freedom in our 
detection scheme as follows. First, we need to introduce the abil- 
ity to measure each mode in an arbitrary phase quadrature at the 
MPHD. This may be achieved by manipulating the phases between 
the copropagating entangled modes. By introducing a unitary 
mode-shaping device such as that explained in references^^'^^, such 
access to individual phases becomes possible. These papers describe 
a process of manipulating both the amplitudes and phases of 
spatial modes in a lossless fashion via an adaptive process utilizing 
a network of deformable mirrors and lenses. 

Second, in order to perform teleportation-based computations, 
we require the ability to feed-forward the measurement results of 
arbitrary cluster nodes to remaining cluster nodes. This may be 
realised by combining a specifically shaped displacement beam 
containing the necessary phase space displacements with the beam 
containing the copropagating modes. By shaping the displacement 
beam to spatially match a specific mode to be displaced, it will inter- 
fere only with the desired mode owing to the orthogonal nature 
of the spatial mode basis. These modifications are feasible with 
existing technologies^^. 

Emulating linear optics networks by mixing copropagating 
spatial modes is a highly efficient method for generating multi- 
mode entanglement. Otherwise arduous and potentially lossy tasks 
such as mode matching during the construction of a linear optics 
network are performed effortlessly and losslessly via software- 
controlled combinations of the spatial modes. We have shown that 
although correlations weaken if more squeezing resources are not 
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Figure 5 | Inseparability for different entangled mode bases. The 

solid black line represents the bound of separability. Dashed lines 
represent theory. All uncertainties for experimental values are ±0.01. 
(a) The blue markers are the averaged measured experimental values 
for N-mode inseparability (right column of Table 1), and the dashed 
blue line joins the theoretical values of inseparability with the same 
two squeezed inputs used in the experiment. All experimental losses 
have been taken into account. The red circles are the measured 
experimental values for N-mode cluster states with theory indicated 
again by the dashed red line. Here the maximum value of each row in 
Table 2 is shown rather than the average value, in order to show that 
cluster states have a much more stringent requirement on squeezing 
levels (N = 5 is clearly separable and not a cluster state), (b) All traces 
have two squeezed inputs and N-2 vacua modes, as in the experiment. 
What changes is the amount of squeezing in the two squeezed inputs, 
assumed here to be symmetric with equal antisqueezing. From the 
top we have: -1 dB (magenta); - 3 dB (cyan); experimental parameters 
(blue); experimental values (blue markers); -6dB (green); and 
-IOdB(red). 



Virtual networks. For even-numbered mode bases {N= 2, 4, 6 and 8), the 

method for creating the virtual network is as follows. The two squeezed modes 
ai and a2 are combined on a half reflecting beamsplitter (HBS). As the output 
of this HBS is an EPR state, we choose to call this the BBS. The BBS outputs are 
symmetrically combined with N- 2 vacua, as in F ig. 1 of the main text. The 
beamsplitters (BSs) are then given — n), where n is the 

number of BSs between the BBS and the BS in question. For N= 4 and N= 6 and 
N= 8, mode output 2 is swapped with mode output N- 1. For N= 8, an additional 
swap of output modes 4 and 5 is made. For odd-numbered mode bases {N= 3, 
5 and 7), the method is the same with the following modifications. The BBS 
has its reflectivity changed to r = 1/2 — 1/{2N) (see, for example, refs 30,31 for 
more details on N= 3). The vacua are mixed using beamsplitters as above, with 
one output arm having one less vacuum input. K phase shifts are applied to all 
BS outputs on the left of the BBS except for the one left output exiting the last 
BS. Mode outputs 1 and N- 1 are swapped, and the network for N= 7 has an 
additional swap between output modes 3 and 4. The homodyne gains gi are 
optimized using a genetic algorithm in order to maximally satisfy the van 
Loock-Furusawa inequalities. These gains gi scale the contributions of the 
quadrature variances and are independent from calculations regarding (ynet- 
Here, optimal homodyne gains are calculated using two measures: minimizing 
the mean of the N-1 inequalities; and minimizing the variance of the set of 
inequalities. A trade-off" between the two measures is needed, and preference 
is given to minimizing the mean of the inequalities. 

Spatial mode bases. The input matrix is defined as follows: 



1 



(5) 



match the detection basis to the input modes. For example, the top row contain- 
ing all ones recovers the standard GM, and the second row recovers the phase- 
flipped FM. By setting (y^et = can label each row of as GjJ^. Formally, 

,(8-N))^in 0N,{8-N))Uin> where 0]v,(8-N) is a zero matrix of 
sizeNby (8-iV). 

The linear optics network for the ideal and symmetric 2 -mode BPR basis is 
simply a HBS: 



■^net ■ 



( 1 


1 


V2 


^/2 


1 


1 


k42 


V2 



(6) 



added, (non-cluster) entangled modes scale here as the number 
of orthogonal modes measurable within the beam. The maximum 
number of measurable modes corresponds directly to the number of 
photodiodes in each pair of the multi-pixel detectors. We have dem- 
onstrated this by measuring iV=2, 3, 4, 5, 6, 7 and 8-mode entan- 
glement within one beam, including up to 4-mode cluster states, 
switching between them in real time. The ability to perform a wide 
range of protocols and optimize networks for asymmetry using just 
one optical setup offers versatility to future networks that will utilize 
entanglement as a resource. 

Methods 

Experimental setup. We use a dual-wavelength continuous-wave Nd:YAG laser 
at 1,064 nm and 532 nm. The OPAs each contain a periodically poled KTP crystal 
in a bow-tie cavity. The squeezed beams are almost identical in purity, with 
squeezing levels of — 6 dB and antisqueezing of 8.5 dB. The beam containing 
the eight spatially orthogonal modes (see main text) is made highly elliptical in 
order to be measured by the MPHD, which has a linear array of eight photodiodes. 
The photodiode array used is a Hamamatsu InGaAs PIN photodiode array 
(G7150), which actually has 16 photodiodes. However, we choose to use only 
eight of these in the present experiment. The filling factor for the array is 90%, 
meaning that 10% of the light does not hit an active surface. The quantum 
efficiency for the photodiodes is 80%. 



From equation (4) we get: 



f .2\ 



1 



1 



72 V2 
1 _ 1 

72 >/2, 



^J_^^/2 ^/2 ^/2 ^/2 0 0 0 0 
Vsto 0 0 0 ^ ^ ^ ^ 



from which we see that indeed Gi and share no part of the detected light. 
We show this ideal BPR basis in Fig. 2 in order to emphasize the spatial 
separation. Note that the factor 1/ Vs has been omitted from the scale in 
Fig. 2 for clarity. 
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The optimized network uses a beam-splitter reflectivity of 48.8%, and produces 
the following output modes: 



.2,opt 



1^2 



2,opt 



2,opt 2J 



1 f 0.699 
" VS t 0.716 

1 f 1.414 
"Vst 0.017 



f A2,opt 



0.716 
-0.699 



1111 
1 -1 -1 -1 



1.414 1.414 1.414 -0.017 -0.017 
0.017 0.017 0.017 1.414 1.414 



-0.017 
1.414 



-0.017^ 
1.414 



A2,opt 
V^2 



Measuring cluster states. The difference between a 2-mode cluster state and a 
2 -mode EPR state is a Fourier transform on one mode. The Fourier transform is 
a rotation of 71 / 2 degrees: 



F = R 



K 71 ^ 

cos(— ) -sin(— ) 
2 2 

sin(-) cos(-) 
2 2 



0 -1 

1 0 



(7) 



Therefore, we get follows that the homodyne measurements 

we perform in the 2-mode cluster basis {pi -^2'p2 "■^ilcluster the 2-mode EPR 
basis {xi-X2,pi +P2}epr equivalent. Therefore, we may perform local Fourier 
transforms so long as we can match the homodyne detection basis for individual 
modes. It is important to note that this convenient basis change will not always be 
possible for different clusters. However, by shaping the local oscillator we may have 
access to arbitrary cluster states within the one beam. This was out of scope for the 
current experiment. 

The criteria for verifying the measurements of the various cluster states are 
given below^^'^^ with the results summarized in Table 2: 



(8) 



(9) 



(10) 



^- = 2 

I:([A(pi-X2)f> + <[A(p2-^l)f><l, 



N = 3 

I:([A(pi-X2)f> + <[A(p2-Xi-X3)]2><l, 
II : ([ A(P2 - ^1 - % ) f > + <[ A(P3 - X2 ) f > < 1 , 



I : <[ A(pi - X2 ) ]^ > + <[ A( P2 - ^1 - ^3 ) ]^ > < 1 ' 
II : ([ A(P2 - ^1 - % ) ]^ > + <[ A(P3 - X2 - X4 ) ]2 > < 1 , 
lll:{[A{p^-X2-h)f} + {[A{p4-h)f)<h 
N = 5 

I : ([A(pi -X2)f} + {[A{p2 -xi-x^)f}<h 
II : ([ A( P2 - ^1 - % ) ]^ > + <[ A(P3 - X2 - X4 ) ]2 > < 1 , 
III : ([ A( P3 - X2 - X4 ) ]2 > + ([ A( P4 - X3 - X5 ) ]2 > < 1 , 
IV : ([ A( P4 - X3 - X5 ) ]2 > + ([ A( P5 - X4 ) ]2 > < 1 . 
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